Tristability in the pendula chain 
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Experiments on a chain of coupled pendula driven periodically at one end demonstrate the ex- 
istence of a novel regime which produces an output frequency at an odd fraction of the driving 
frequency. The new stationary state is then obtained on numerical simulations and modeled with 
an analytical solution of the continuous sine-Gordon equation that resembles a kink-like motion back 
and forth in the restricted geometry of the chain. This solution differs from the expressions used 
to understand nonlinear bistability where the synchronization constraint was the basic assumption. 
As a result the short pendula chain is shown to possess tristable stationary states and to act as a 
frequency divider. 
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I. INTRODUCTION. 

The sine-Gordon model and its discrete analogue, the 
Frenkel-Kontorova chain, are among the most prominent 
equations of nonlinear physics, and have attracted inter- 
est of people working in quite different fields, see e.g. 
books [I], [2, HI • in particularly, topological (kinks) and 
nontopological (breathers) solutions of the sine-Gordon 
equation describe the dynamics of nonlinear excitations 
in various spatially modulated systems, e.g. dislocations 
in crystals [J], magnetic and ferroelectric domain walls 
motion [f|, vortices in arrays of Josephson junctions Q, 
etc. At the same time the model has a simple experi- 
mental counterpart, namely the chain of linearly coupled 
pendula 0] which offers an interesting opportunity to eas- 
ily visualize all the main nonlinear characteristics of the 
sine-Gordon system. Then this simple laboratory tool al- 
lows to observe novel effects @, H, B Gil which may then 
apply in completely different physical situations. 

As a matter of fact, a recent experimental discov- 
ery of supratransmission effect in the pendula chain (Tl| 
has led to the study of similar phenomena in optical 
Brag g g ratings [l^ . Josephson Junction transmission 
line |l3l | waveguide arrays [lj, By this approach, 
many simil ar p henomena observed in the same systems 
[lrl Uj\, [l8l. Il9l. [20| have been identified as effects of non- 
linear bistability. Moreover it allowed us to predict the 
existence of bistable magnetization profiles in thin mag- 
netic films [2l[ and to suggest ultrasensitive detectors 
(or digital amplifiers) in optical waveguides [22| quan- 
tum Hall bilayers [23j and Josephson Junction parallel 
arrays [24| . 

The bistability property can be simply formulated by 
saying that pfjj a given periodic boundary driving may 
produce two completely different stationary states: one 
which tends to the linear evanescent profile at vanishing 
amplitude, the other one which can be qualitatively un- 
derstood as a portion of the stationary breather-like so- 
lution which exists only if the system size is comparable 



with the characteristic length of the fundamental (con- 
tinuous) breather solution. This is a main difference with 
most of the earlier studies on sine-Gordon model where 
the semi-infinite chain has been examined (sec e.g. Rcf. 




FIG. 1: Pictures of three stationary states of the pendula 
chain obtained for one single given driving amplitude and 
frequency. The upper graph corresponds to the lowest energy 
state and the lower one describes the largest energy regime. 
An approximate energy hierarchy is 1:10:100. 

We report here the discovery of a third stationary state 
which can be qualitatively understood as the motion back 
and forth of a kink-likc structure in the short pendula 



2 



chain. The new stationary regime apears to be com- 
pletely different from the two cases considered earlier e.g. 
in [HI, [13, [Hill HI - As a matter of fact, such a dynam- 
ics creates a new frequency in the system and furnishes 
a tool to divide the input frequency by odd fractions (we 
shall illustrate chain end oscillations with frequency f2/3 
or f2/5, where fl is the driver frequency). The value 
of the odd divider depends both on the input frequency 
range and on the length of the chain. Let us recall that 
the previously discovered two regimes are synchronized 
to the driver (same input and output frequencies). 



II. MODEL EQUATIONS. 

The dynamics of the chain of N pendula is naturally 
described by the Frenkel-Kontorova model [l[ 



+i 



2u n ) + sinu„ = 0, (1) 



where overdot means derivation with respect to time. 
The variable u n is the angular deviation of the n th pen- 
dulum, ujq is the cigenfrequency of a single pendulum 
and a is proportional to the linear torsion constant of 
the spring (for our experimental chain luq = 15.1 Hz, 
(7 = 32.4 Hz). The damping coefficient 5 is phenomcno- 
logical, it has been evaluated in the experiments as ap- 
proximately 8 = O.Olwo- This is the value actually used 
in the numerical simulations. The applied periodic driv- 
ing is here modeled by the boundary conditions 



uo(t) = b cos(fit), un + i = un 



(2) 



which model a forced end in n = and a free other end 
in n = N. 

It is worth insisting on the fact that the chain is sub- 
mitted to a prescribed boundary value (the datum of 
uo(t)), not to a given force acting on the first particle. 
In the experimental setup, the motion of the virtual pen- 
dulum uq is the driving engine motion which has indeed 
a prescribed motion obtained through feedback control. 
This has an important fundamental consequence: the en- 
ergy absorbed by damping is continuoulsly compensated 
by the action of the driving and therefore the motion of 
the chain is quite similar to that of an undamped device 
(for which, in a stationary regime, the driving would not 
give energy to the chain). 

The experiment consists thus in driving the short chain 
pictured in Fig[I]with a frequency in the forbidden band 
gap (CI < luo), which actually does not excite linear 
modes. Without external perturbation the system locks 
to a periodic solution with low output amplitude wjv(i). 
Depending on the value of an external kick one makes the 
system bifurcate to two different stationary states and we 
thus obtain a tristable behavior (with approximate hier- 
archy of energies 1:10:100). The reader will find on the 
web page (27j a movie of the experiment where the sys- 
tem is first set in the high-energy new stationary stable 
regime and then put down successively to the two others 
stable states by taking energy off. 
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FIG. 2: Analytic input-output amplitude dependences for dif- 
ferent oscillation frequencies of stationary states given by for- 
mulas ((5]l where continuous line represents solution (I), dot- 
ted line solution (II) and dashed line solution (III), plotted for 
two frequencies, namely 0.9 ujq and 0.3 ujq as indicated. The 
points 1, 2 and 3 corresponds to the stable regimes with a 
single driving amplitude \uo(t)\ max = 0.5 rad. The points 1 
and 2 represent the situations when the whole chain oscillates 
with the driving amplitude 0.9 u>o but with different output 
amplitudes. The point 3 corresponds to the driving frequency 
0.3 wo and describes kink motion forth and back. As the ex- 
periments and numerical simulations show (and this is a main 
finding of the paper), the latter regime can also be reached 
with a driving frequency 0.9 loq three times larger than the 
one actually used. 



To develop an analytical description of the process, let 
us consider the continuous approximation of eq.JT]) by 
substitutions t — > u>ot, n = loqx/ 'a. Neglecting dissipation 
we obtain the sine-Gordon equation 



e [o,L] 



Utt 



s'mu = 0, 



(3) 



where L = Na/u>o. The mixed Dirichlet and Neumann 
Boundary Conditions u(0, t) = b sin(rtt) (driven bound- 
ary), u x (L,i) = (free end boundary) allows to seek the 
following periodic stationary solutions p5| 



u(x, t) = 4 arctan 



X(x)T(t) 



(4) 



where one has three choices (cn, sn and dn are the stan- 
dard Jacobi elliptic functions) 

(I) X = cn(f3(x - L),n), T^cn(ujt,v), 

(II) X = dn(0(x -L),fj), T = sn(ut, u), (5) 

(III) X = da.(J3(x — L) + K(ji), fx), T=m(ut,v). 

Here K(/i) stands for a complete elliptic integral of the 
first kind of modulus fi. These families of solutions are 
parametrized by the two free constants ui and v £ [0, 1], 
then for solutions of type (I) the remaining parameters 
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are given by 

l 4 2/i 2\ 2 2 



2r = 1 - u 2 + 2wV + y/(l -co 2 ) 2 + 4lu 2 is 2 , 
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r b + r 2 ' 

while in both cases (II) and (III) they read 

7 4 2 2 2 

= 0;^, s = — uj v , 



(6) 



2r = 1 - oj 2 (1 + v 2 ) + - oj 2 (1 + v 2 )} 2 - Auj^v 2 



(7) 



Note that r should be real valued and positive which may 
restrict the allowed values of u. 



III. TRI-STABILITY AND FREQUENCY 
DIVISION. 



Since the experiments (confirmed by numerical simu- 
lations later on) show that the frequency 0/3 can also 
be excited, we assume that the period of the time depen- 
dent part T(t) of the stationary solutions §5§ coincide 
with an odd integer fractions of the driving frequency f2. 
Recalling that the period of T(t) is 4K(u)/cj we require 
thus 



lu = 2VLK(v) / {mir) 



(8) 



where m is an odd integer. For a given value of the pa- 
rameter v g [0,1], the above relation fixes the second 
parameter to in terms of the driving frequency Q. There- 
fore fixing Q (driver frequency) and varying v one can 
plot the output amplitude u(N, t) in terms of the input 
u(0,t) from the analytic expressions We display this 
dependence for SI = 0.9 (in units of u>o) as a the full 
line in figH] where different colors report to different so- 
lutions. We also plot (dashed line) the output amplitude 
for a driving frequency Q = 0.3. Therefore, to the given 
driver amplitude max t |u(0, i)| = 0.5 may correspond two 
stable synchronized states (points 1 and 2 on the graph) 
having the driver frequency 0.9 and one more stable state 
with frequency 0.3. 

It is then a simple matter to check that the stationary 
state related to point 3 of the plot of fig[2] corresponds 
effectively to our numerical simulations, and hence to 
the experiments of figfl] It is done in figO where the 
last plot shows the result of a numerical simulation (full 
line) compared with the analytic solution (dahsed line) 
related to point 3 of figH] Wa have also plotted the time 
evolution of the total energy of each pendulum given by 
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= 2^n + ~T [( u n+l - Un) 2 + («n-l 



+ UJ 2 (1 - COSU n ). 



(9) 
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FIG. 3: Numerical simulations on Frenkel-Kontorova model 
(TTJ with a damping constant S = 0.01 ujq and 8 pendula. The 
time evolution of pendula energy and input-output oscilla- 
tions are displayed corresponding to the points 1), 2) and 
3) in Fig. 2. The driving amplitude is |uo(i)|maa; = 0.5rad 
and its frequency Q = 0.9 u>o for all 3 cases. This results in 
the same output frequency oscillations Q in graphs 1) and 2) 
but f2/3 output oscillations in graph 3). Dashed lines dis- 
play analytical curves obtained from ([5]), while dotted-dashed 
and solid lines represent time evolution of input and output 
oscillations, respectively. 



The numerical simulations of the process are done by 
applying to the model ([1]) the boundary conditions @ 
with b = 0.5 and u = 0.9 luq, together with an initial 
condition where a few pendula at the end of the chain 
are given large initial amplitude. For instance, to reach 
the new stationary state (3) of figGS the chosen initial 
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amplitude is 2w, while for the value 7r, the system locks 
to the state (2). 

It is worth noting that both experiments (as those dis- 
played in [27| ) and numerical simulations contain intrin- 
sic damping. Still the analytic solutions of the continuous 
undamped sine-Gordon model fit strikingly well numeri- 
cal simulations of the discrete damped Frenkel-Kontorova 
model ([T]) . This is a general property of such short length 
driven systems to lock on fundamental solutions of the 
undamped limit, as previously displayed in p3l. [241. |28| . 
The main fact is that, without damping, the chain in a 
stationary regime does not absorb energy and the bound- 
ary value does not transfer any power to the chain. With 
damping, the driving boundary does transmit power to 
the chain in such a way as to compensate exactly the 
losses. The point is that the system is submitted to a 
prescribed boundary value which adapts to the amount of 
lost power and keeps amplitude and frequency constants 
(in the experiments presented in fig[TJ the engine has a 
feedback driving mechanism that controls the amplitude 
and frequency). Last, the proof that the analytical solu- 
tions constitute an attractor for the damped system is, 
as far as we know, an open question. 
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Thus we have actually demonstrated the possibility of 
conceiving a frequency divider with which the driving 
frequency can be divided by 3, 5, 7, depending on the 
chain length. For example the figH] shows numerical sim- 
ulations on a chain of 12 pcndula with a resulting fre- 
quency division by 5 at a driving frequency 0.9 loq and 
by 3 at a driving frequency 0.5 ujq. In such a case we 
have obtained that the frequency is divided by 5 if the 
driving frequency f2 is in the range 0.88w < ft < 0.92w , 
while the same chain can divide the frequency by 3 when 
0.41w < ft < 0.6w . 



IV. CONCLUSION. 

For a single monochromatic driving (fixed ampli- 
tude and frequency), we have demonstrated experimen- 
tally and numerically the existence three states which 
have been given analytic expressions (in the continuous 
limit): the first one is the quasi- linear solution (actu- 
ally a breather- like tail), the second one resembles half 
a breather, both of them oscillating with the driving fre- 
quency f2, and which were already known as the build- 
ing blocks of nonlinear bistability. The discovered third 
state resembles a kink moving back and forth with the 
frequency These sates have been given explicit ana- 
lytical expressions in the continuous limit: the first state 
is described by the solution of type (III) in (JSJ) while 
the solution of type (I) describes altogether the "half- 
breather" with frequency Q and the "oscillating kink" 
with frequency f2/3. 

The process of frequency division is thus induced by 
the motion back and forth of a kink-like structure inside 
the chain. It is possible to extend these studies to other 
realistic physical systems governed by the sine-Gordon 
equation. We expect such a new stationary stable regime 
to be interesting for applications where one is interested 
in producing an odd fraction of the driving frequency. 
The fraction number depends both on the length of the 
chain and on the input frequency range. 

Last but not least, many other well known nonlinear 
systems exhibit nonlinear bistable behavior, as e.g. the 
nonlinear Schrodingcr equation or the coupled mode sys- 
tem in Bragg media, and this discovery is very likely to 
apply also there. 
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FIG. 4: Time evolution of pendula energy and input-output 
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the input frequency when the input frequency is fl = 0.9 wo 
(upper graph). In the lower graph the driving frequency is 
Q = 0.5 u>o and one has frequency division by factor 3. 
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